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Abstract 



We have studied two-body charmless decays of the B meson into the final 
states p°p°, K*°p°, K*°K*°, K*°K*°, K*+p°, K* + K*°, and K*+K*~ using 
only decay modes with charged daughter particles. Using 9.7 million BB pairs 
collected with the CLEO detector, we place 90% confidence level upper limits 
on the branching fractions, (0.46 — 7.0) x 10~ 5 , depending on final state and 
polarization. 
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In the Standard Model, CP violation is introduced by the complex phase in the Cabibbo- 
Kobayashi-Maskawa quark-mixing matrix. The experimental study of CKM phases will 
probe the Standard Model description of CP violation. This may provide a window to 
new physics. In particular, it has been suggested that we may construct a relationship 
between charmless B — > VV decays that may lead to the extraction of the angle a. Earlier 
observations of rare charmless decay modes at CLEO include B — > K-k, titt, t]K, pit, rj'K, 7]K* 
and urn [[§. It is natural to extend our search toward other rare charmless B decays. 

In this letter, we present results of searches for B meson decays into the vector mesons 
p°, K*° and K* + . The decays are dominated by the b —>■ u tree-level and b —> dg penguin 
processes, though other mechanisms may also contribute ||. 

The data used in this analysis were collected by the CLEO detector ||] at the Cornell 
Electron Storage Ring (CESR). The data consist of an integrated luminosity of 9.1 fb _1 at 
the T(4S) resonance, corresponding to 9.7 x 10 6 BB events. To determine backgrounds due 
to non-resonant e + e~ — > qq process, we also collected 4.6 fb _1 of continuum data at energies 
just below the T(AS) resonance. 

The CLEO detector has 67 tracking layers and a Csl electromagnetic calorimeter that 
provides efficient n° reconstruction, all operating within a 1.5T superconducting solenoid. 
The central tracking system, consisting of an inner 6-layer straw tube precision tracker, a 
10-layer vertex drift chamber, and a 51-layer main drift chamber, provides a measurement 
of momenta of charged particles and the vertex position of decaying K$- It also measures 
the specific ionization loss, dE/dx, which is used for particle identification. The precision 
tracker was replaced by a silicon vertex detector for the latter 65% of data taking. Muons 
are identified using proportional counters placed at various depths in the steel return yoke 
of the magnet. 

B candidates are selected by straightforward criteria based on energy-momentum conser- 
vation and event shape. Simulations of the signal and backgrounds are used to refine these 
criteria and to determine their effectiveness. 

The B — > VV decays are reconstructed through the decay channels B — ► p°p°, B° - 
K*°p°, B° -»■ K*°K*°, B° -»■ K*°K*°, B + -> K* + p°, B + -> K* + K*°, and B° -> K*+K*~. 
We form p° candidates from tt + 7t~ pairs with an invariant mass within 150 MeV/c 2 of the 
nominal p° mass. K*° / K*° / K*^ candidates are selected from K^ir^ / 'K^^ pairs within 50 
MeV/c 2 of the nominal K* mass. 

Charged tracks are selected by requiring them to pass quality criteria and must be consis- 
tent with production from the primary interaction point (except for pions from Kg decays). 
The measured specific ionization (dE/dx) of charged kaon and pion candidates is required 
to be within 3.0a (standard deviation) of their most probable values. We reject electrons 
based on dE/dx and the ratio of the track momentum to the associated shower energy in 
the Csl calorimeter. We reject muons by requiring that the tracks not penetrate the steel 
absorber past a depth of 3 nuclear interaction lengths. The Kg is selected by requiring a 
decay vertex displaced from the primary interaction point and an invariant mass within 10 
MeV/c 2 of the K s mass. 

Fully reconstructed B mesons are selected on the basis of the beam-constrained mass of 

the candidate, Mb = y E% eam — Preconstructedi anc ^ ^ e difference between the reconstructed 
and beam energies, AE = E reconstructed — E beam . AE is sensitive to missing or extra particles 
in the B candidate, as well as incorrect assignment of particle masses. For the fully recon- 
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structed B meson decays in this analysis, the Mb distribution peaks at 5.28 GeV/c 2 with 
a resolution ranging between 2.2-2.6 MeV/c 2 , and AE peaks at zero GeV with a resolution 
ranging from 16 MeV to 27 MeV. Candidates are accepted for further analysis if AE and 
Mb are within a signal region ±2cr around the central signal values for all channels (except 
K*~K* + where a larger region is used since this involves two K s 's and is therefore relatively 
clean). 

The backgrounds consist primarily of continuum events from e + e~ — > qq (q = u, d, s, c) 
with a 10-15% contribution from B decays, and are estimated from a combination of off- 
resonance data and b — > c Monte Carlo. Event-shape variables can be used to discriminate 
against the jet-like continuum events since B mesons are produced nearly at rest. Accord- 
ingly, we select only events with R 2 < 0.5, where R 2 is the ratio of the second to zeroth 
Fox- Wolfram moments of the event H . In continuum events, momentum conservation aligns 
the thrust axis of the B candidate with that of the rest of the event while they are almost 
uncorrelated in BB events. This allows additional suppression of continuum by restricting 
| cos #tt |, the angle between the two axes. We require |cos^| < 0.7 for all decay modes, 
except for K* + K*~ , where we use | cos6 tt \ < 0.9. 

The four selection criteria discussed above, on M B , AE, R 2 and cos9 tt , determine the 
signal efficiency (e) for each mode. We measure this efficiency using Monte Carlo simulation 
for each of the 3 possible helicity states of the decay products: 00, -1-1 and +1+1. Our 
study indicates that the 00 helicity has slightly lower efficiency than the 11 helicities, since 
it results in more low momentum charged pion and kaon tracks from the B decay chain, 
for which the detector has a lower acceptance. In addition, the 00 state will tend to align 
the vector decay products leading to a higher average R 2 , also decreasing the efficiency. 
We give separate results assuming the signal is 100% 00 helicity or 100% 11 helicity. For 
any assumed helicity distribution of signal events in the data sample, upper limits can be 
obtained by linear interpolation. 

We find significant double counting of events in the K*°p° channel, caused in most cases 
by the Kj tt ambiguity in the K*° — ► K + + it~ sub-decay. In the final results we count only 
one entry for each event. We also consider the possibility of cross-feed between different 
channels of B — > VV decays. Neglecting the contribution from the forbidden decay mode 
B — > K*°K*° (AS = 2), the cross-feed effect is small even if we use the 90% upper limits 
to evaluate the cross-feed contribution to the yields. We do not correct for this contribution 
when extracting the upper limits. 

There are several sources of systematic error. A substantial contribution comes from the 
uncertainty in track efficiency, which is 1.5% per charged track. For B decay modes with 
K*^, there is an additional 5% uncertainty due to the Kg vertex requirement. In addition, 
we estimate 1% per charged track uncertainty due to the dE/dx requirement. Additional 
systematic errors include 7% uncertainty from the thrust criterion and 3% from the AE 
and Mb requirements. Uncertainties due to Monte Carlo statistics range from 2% to 6%, 
depending on B decay mode. 

The results of this analysis are summarized in Table 1; we see no statistically compelling 
signal in any individual decay channel. To calculate 90% confidence level (C.L.) upper limits 
on the number of signal events in each channel, we use the Poisson likelihood of a hypothesis 
for the average number of signal events ns given n b s events detected and a background of 
n B = n 6 _ c + n ff. 
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TABLE I. The 90% C.L. upper limits for the B — > VV decay modes (£>cleo) are shown in units 
of 10~ 6 , along with the corresponding theoretical predictions (^theory) ||- n obs is the number of 
observed events, n tt is the off-resonance background (normalized), n(,_> c is the BB background 
estimate (from Monte Carlo), and n u x. is the corresponding upper limit including systematic error 
and background statistics. The reconstruction efficiency (£) is also shown along with the systematic 
error (5e). We assume equal branching fractions for T(4S) — > B°B° and B + B~ . 
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C(n s ,n B ,n obs ) = e^ ns+nB \n s + n B ) n °»° / n obs \ 

To include the effect of the systematic error in the acceptance and the limited statistics of the 
background samples, we convolute the appropriate distributions (Gaussian for the systematic 
error, and Poissons for the h — > c and continuum backgrounds) with the likelihood function, 
C, to obtain a modified likelihood function, C* . The 90% C.L. upper limit on the yield is 
obtained by integrating £*, i.e. by solving for n u j, in: 




C*(n s ,n B ,n obs )dn s = 0.90 



The upper limits on the branching ratios are then calculated from the formula, 

B(B -> VV) = — 

n B B x e x 11b 

where n BB is the number of BB meson pairs in the data sample, and Y\ B is the product over 
all the relevant branching fractions of the vector meson decay chain. 

To summarize, we set 90% C.L. upper limits on branching fractions of seven B —>■ VV 
charmless decay modes. Theoretical predictions for the branching fractions of these modes 
tend to be near 10~ 6 . Thus our results are consistent with theoretical calculations based on 
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the Standard Model. In order to challenge these predictions data samples of the order of 10 8 
BB mesons would be required. 
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